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ABSTRACT 

We examine the accretion and merger histories of central and satelhte galaxies 
in a smoothed particle hydrodynamics (SPH) cosmological simulation that resolves 
galaxies down to 7 x IO^M^q. Most friends-of- friends haloes in the simulation have a 
distinct central galaxy, typically two to five times more massive than the most massive 
satellite. As expected, satellites have systematically higher assembly redshifts than 
central galaxies of the same baryonic mass, and satellites in more massive haloes form 
earlier. However, contrary to the simplest expectations, satellite galaxies continue to 
accrete gas and convert it to stars; the gas accretion declines steadily over a period 
of 0.5 — 1 Gyr after the satellite halo merges with a larger parent halo. Satellites 
in a cluster mass halo eventually begin to lose baryonic mass. Since z = 1, 27% of 
central galaxies (above 3 x lO^^M©) and 22% of present-day satellite galaxies have 
merged with a smaller system above a 1:4 mass ratio; about half of the satellite 
mergers occurred after the galaxy became a satellite and half before. In effect, satellite 
galaxies can remain "central" objects of halo substructures, with continuing accretion 
and mergers, making the transition in assembly histories and physical properties a 
gradual one. Implementing such a gradual transformation in semi-analytic models 
would improve their agreement with observed colour distributions of satellite galaxies 
in groups and with the observed colour dependence of galaxy clustering. 

Key words: galaxies: evolution — galaxies: formation — models: semi-analytic — 
models: numerical 



1 INTRODUCTION 

In the standard theoretical description of galaxy formation, 
galaxies grow initially by the condensation of gas at the 
centr e s of dark matter potent ial wells (e.g. IWhite fc ReesI 
1 19781 : iFall fc Efstathiou I Il98(il ). When a large dark mat- 
ter halo accretes a smaller halo, the galaxy hosted by the 
smaller halo becomes a satellite system, which may even- 
tually merge with the central object after sinking by dy- 
namical friction. Even if the satellite retains some of its 
original dark matter, it is likely to move through the large 
parent halo with a random velocity that exceeds its own 
escape speed, making it difficult for the satellite to ac- 



crete gas from the halo or to merge with other satef 
lites. In semi-analytic models of galaxy formation (e.g 
White fc Freii^liggil: iKauffmann. White fc Guideronilll993l: 



Cole et all Il994l: iPrimack et all Il998l : ICroton et all l2006l : 
Bower et al. 20061 ) . these expected differences between cen- 
tral and satellite galaxies are usually encoded by simple 
rules, with some variations from one model to another. 
When two haloes merge, the central galaxy of the more mas- 
sive progenitor becomes the central galaxy of the new, larger 
halo. Cooling gas from the halo is assigned mostly or entirely 
to the central galaxy. Satellites accrete little or no fresh gas, 
and their pre-existing gas may be removed by ram pres- 
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sure stripping when they enter the larger halo. Most or all 
mergers are assumed to involve satellites joining the central 
galaxy rather than merging with each other. 

In this paper, we examine the growth of central and 
satellite galaxies in smoothed particle h ydrodynamics (SPH) 
cosm ol ogical simula t ions IIKatz. He rnguist. & Weinberg 



19921 : lEvrard et all 1 1994 
19961 ). In such simulations, differences between central 



Katz. Weinberg fc Hernguistl 



and satellite galaxies are not imposed a priori, but central 
galaxies nevertheless emerge as a distinct class for the 
physical reasons discussed above. The galaxy closest to the 
halo centre of mass is usua lly older and more massive than 
satell ites in the same halo l|Berlind et al ] |2003l : IZheng et al.l 
I2OO5I ). Here we investigate the assembly histories of present 
day centrals and satellites as a function of galaxy and 
halo mass, and we track the growth of galaxies to see the 
extent to which satellite systems experience continuing gas 
accretion or mergers with other satellites. 

Our goals are to understand the physics of galaxy as- 
sembly in the simulation and to inform semi-analytic mod- 
els of galaxy formation and other efforts to interpret ob- 
served galaxy clustering and the environmental dependence 
of galaxy properties. The central-satellite distinction is a 
key element in modeling the origin of bimodality in the 
galaxy population, in which galaxies with old stellar pop- 
ulations and little ongoi ng star formation form a distinct, 
passive, "red sequence" (fstrateva et al.l 200 ll : iBlanton et al.l 
I2OOII : iBell et al.ll2003l : iFaber et al.l 120071 ). Strangulation of 
gas supplies in satellites is the likely route by which many 
galaxies join the passive sequence, though additional mech- 
anisms are required to explain the red colours and low 
star formation rates of the most massive galax:ies (see e.g. 
Croton et al. I I2OO6I : iDekel fc Birnboim|[200al : ICattaneo et al.l 



20071 ). Central-satellite separation is also a key ingredient 



in "halo oc cupation distribution" (HOD) models o f galaxy 
clustering (|Kravtsov et all 12004 IZheng et all |2005| ) and in 
recent efforts to i nterp r et properties of gala-xy groups (e.g. 
Woinmann ct al.' '2006'; 'Cilbank et alj I2OO8I : iHansen et all 
2007; Li ct al. 2007; Coil ct al. 2008). 

The simulation th at we analyse here is a lso the main 
simulation analysed bv lKeres et aU l|2005l ) and lMaller et al.l 
l|2006l ). who examine gas accretion and merger rates but 
do not dist inguish between central and satellite galaxies. 
iKeres et al.l (fioOSi ) emphasise the distinction between be- 
tween "cold" and "hot" modes of gas accretion. In haloes 
of total mass M < 3 x W^^AIq, most gas that accretes 
onto galaxies does so without ever heating close to the 
halo virial temperature. In higher mass haloes, most gas 
heats to the virial temperature before settling into galax- 
ies. Overall, about half of the baryonic mass in the sim- 
ulation's z = galaxy population was originally accreted 
in "cold mode" and half in "hot mode", with the former 
dominating at high redshift and in low mass galaxies today. 
Here we will investigate both total and cold-only accretion 
rates of central and satellite syste ms. Further di s cussions 
of th e cold/hot dichotomy appear in 'Binnov' (1977"), BinnevI 
1(2004 ) ■ Katz (1992), Kay ct al. (2000), Fardal ct al. (20oC 
iKatz et all (|2003|). iBirnbqim fc Dekell ||2003|). ICroton etal] 



l 2006t ). lDekel fc BirnboimI l|2006l ) and lKeres et al.l (|2008l ). 

We describe our simulation and methods for identifying 
haloes and galaxies in §2, before proceeding to our main re- 
sults in §3. There we start with a qualitative overview of ac- 



cretion and merger properties of central and satellite galax- 
ies, then quantify their relative accretion rates and merger 
rates and their dependences on halo mass. In §4 we sum- 
marise our results and discuss their implications, with par- 
ticular attention to recent studies of galaxy bimodality and 
the red population in groups. 



2 SIMULATION AND NUMERICAL 
METHODS 

2.1 Simulation 

Our simulation i s performed using a parallel implementa- 
tion of TreeSPH (jP ave ct al."l997|; iHemauist fc Katz|[l989l : 
iKatz. Weinberg fc He rnguist 1996i) : which combines a tree 
algorithm for gravitational calculations with smoothed par- 
ticle hydrodynamics (SPH). There are three kinds of parti- 
cles in our simulation: dark matter, stars and gas. The coUi- 
sionless particles (dark matter and stars) are influenced only 
by gravity, whereas the gas particles are influenced by pres- 
sure gradients and shocks in addition to gravity. Gas parti- 
cles experience adiabatic heating and cooling, shock heating, 
inverse Compton cooling off the CMB and radiative cooling 
via free-free emission, coUisional ionisation and recombina- 
tion and coUisionally excited line cooling. Further details 
of the code can be found in iKatz. Weinberg fc Hernguistl 
1 19961 ). 

We model a comoving 22.222h~^ Mpc periodic box with 
128'^ gas particles and 128'' dark matter particles. The dark 
matter particle mass is 7.9 x 10* Mq , and the SPH particle 
mass is 1.05 x 10* M©. The gravitational force softening is a 
comoving 5h~^ kpc cubic spline, which is roughly equivalent 
to a Plummer force softening of 3.5h~^ kpc. 

We include star fo r mation as described by 
iKatz. Weinberg fc Hernguistl (|l996l ). Overdense, Jeans 
unstable gas that is part of a converging flow with density 
Pgas > O.lrriHcm"* and temperature T < 30, OOOK is 
converted to stars on a timescale set by the dynamical 
timescale or the cooling time, whichever is longer. Our 
prescription for star formation leads to a relation with 
g as surface dens i ty sim ilar to a Schmidt law. We adopt 
a iMiller fc Scald (Il979l ) initial mass function. We include 
supernova feedback with 7.35 x 10~^ supernovae per solar 
mass, and each supernova deposits 10^^ ergs of energy 
into the surrounding medium. Since the surrounding gas is 
dense, this energy is usually radiated away before it can 
drive outflows or suppress subsequent star formation. 

We adopt a ACDM cosmology (inflationary cold dark 
matter with a cosmological constant) with flm~0.4, S7a=0.6, 
/i=0.65, Qb = 0.02h~^, spectral index n=0.93, and a8=0.8. 
These values are reasonably clos e to the estimates from 
the cosmic microwave background (|Dunklev et al.ll2008l ) and 
large scale structure (|Tegmark et al.l 20061 ) . though these 
favour somewhat lower ilm (~ 0.25) and higher h (~ 0.7). 
We do not expect modest differences in cosmological param- 
eters to influence our conclusions. 

Hydrodynamic cosmological simulations produce dense 
groups of baryons with sizes and masses c omparable t o 
the luminous regions of observed galaxies l|Katz I Il992l ). 
We identify galaxies using the Spline Kernel In tcrpolative 
DENMAX (SKID) algorithm (|Gelb fc Bertsching cr 199|; 
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iKatz. Weinberg fc Hernauistlll99d ). which identifies gravi- 
tationally bound particles associated with a common density 
maximum. We apply this algorithm to all baryonic particles 
with T < 3 X lO'^K and Pgas/pg^s > 10^ and we refer to 
the groups of stars and cold gas thus identified as galaxies. 
The simulated galaxy population becomes substantially in- 
complete belo w our resolution threshold of 64 SPH particles 
l|Murali et al.L 2002). which corresponds to a baryonic mass 
of 6.8 X 10^ Mq. We ignore lower mass galaxies in our anal- 
ysis. There are 1120 galaxies above the resolution threshold 
in our simulation. 

We identify da rk matter haloes u sing a FOF (friends-of- 
friends) algorithm (|Davis et al.ll 19851 ). The algorithm selects 
groups of particles in which each particle has at least one 
neighbour within a linking length. We adopt a linking length 
I = 0.2n~^/^, where n is the mean dark matter particle 
density, which leads to haloes with typical mean overdensity 
pI'P K, 200. By chance, the simulation contains one halo 
that is anomalously massive (Mh = 3.4 X IO^Mq) for the 
simulation volume, which has some impact on our statistics 
for satellite galaxies when these are not broken down into 
halo mass bins. 

We use 227 output epochs between z = 9 and z = 0, 
with a typical spacing between output epochs of less than 
150 Myr. At each epoch, we identify haloes and galaxies 
above the resolution threshold. 



2.2 Central and satellite galaxies 

Figure[T]shows a selection of six FOF haloes and their galaxy 
populations at z = 0. Each panel is centred on the most 
bound dark matter particle in the halo. The centre of mass 
is shown as a red filled circle in the figure. The galaxies are 
shown as open circles with areas proportional to their (bary- 
onic) mass. A large majority of haloes have a clear central 
galaxy that is more massive than other galaxies in the halo 
and located at the bottom of the dark matter potential well. 
However, in a small number of haloes there are two galax- 
ies of comparable mass, each with its own satellites. Some of 
the systems without a clear central galaxy consist of distinct 
sub-groups that are in the process of undergoing a merger, 
as in the lower left panel of Figure [T] 

We hereafter define the central galaxy of each halo to 
be the galaxy with the largest baryonic mass. Inspection of 
many haloes shows that this identification is always reason- 
able when the halo is regular, and as reasonable as any other 
choice for the small fraction of irregular haloes. 

Figure [5] shows the masses of galaxies and haloes in the 
simulation. The left hand panel plots galaxy mass versus 
halo mass, with central galaxies shown as filled squares and 
satellite galaxies shown as open circles. The right hand panel 
plots the distribution of Mi /Mo, where Mo is the mass of 
the central galaxy and M\ is the mass of the most massive 
satellite. We omit haloes that have no satellites above the 
resolution threshold. The distribution of satellite mass ratios 
is broad, and more than two-thirds of haloes have 0.2 < 
Ml /Mo < 0.5. Only a small fraction of haloes have two 
galaxies of comparable mass (Afi/Mo > 0.8). 
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Figure 1. A selection of friends-of-friends haloes and their galaxy 
populations. The small blue dots are dark matter particles, the 
galaxies are shown as open black circles with their areas propor- 
tional to their masses. Each panel is centred on the most bound 
dark matter particle shown as a blue cross, and the centre of mass 
of the dark matter particles is shown as a red filled circle. 
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Figure 2. 

(Left) Masses of central galaxies (squares) and satellite galaxies 
(circles) plotted against halo mass. (Right) A histogram of the 
ratio of the mass Mi of the most massive satellite galaxy to the 
mass Mo of the central galaxy. 
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Figure 3. Evolution of galaxies with time. Cold accretion is shown in blue and hot accretion in red. M/^ (top left corner) is the halo 
mass. The black crosses show the subsequent evolution of the mass of the galaxy after its last accretion event. For the satellite galaxies, 
the magenta coloured vertical line marks the epoch at which it became a satellite. The first column shows central galaxies, while columns 
2, 3 and 4 show satellite galaxies in decreasing order of halo mass. Panels (a) and (b) show galaxies with masses of ~ 3 X IO^^Mq while 
panels (c) and (d) show galaxies with masses of ~ 10^^ Mq. 
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Figure [Sjcontinued) Panels (e) and {/) show galaxies with masses of ~ 5 X IO^^A/q while panels (3) and (K) show galaxies with masses 
of ~ 2 X IQlf Mq. 



© 0000 RAS, MNRAS 000, 000-000 




0.5 



3 



-0.5 



-1 



I 1 1 1 1 1 > 1 1 1 1 1 1 1 1 1 1 

- t=ll.UC>r 


1 1 1 1 1 1 1 1 1 
















t .... 1 .... 1 .... 1 . 





-1 -0.5 0.5 


1 1.. 



0.5 



-0,5 



K = D 

L-K3.^ Cjr 




Sc/h I Mpc 



-0,5 Q 0.5 



Figure 4. The evolution of a 3 X IO^'^A/q halo with redshift, with positions marked in comoving Mpc. Dark matter particles are 
shown as black dots. The particles in the central galaxy at 2: = are shown in yellow. The red galaxy is the most massive satellite galaxy, 
which becomes a satellite at 2 = 0.514, while the green galaxy becomes a satellite at z = 0.736. The baryonic particles of other z = 
satellites are shown in blue. The inset panels in the right hand corners show a magnification of the red and green satellites, large enough 
to encompass all of their member particles. While the central galaxy has a large total extent at 2 = 0, its stellar half-mass radius is only 
3.9 h-^ kpc. 
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3 ACCRETION AND MERGER PROPERTIES: 
CENTRAL AND SATELLITE GALAXIES 

FigureOshows the assembly history of a representative sam- 
ple of central and satellite galaxies. At each of our 227 output 
epochs from the simulation, we use SKID to identify galax- 
ies with at least 8 SPH particles. For galaxies resolved in 
two consecutive output epochs, we identify gas particles that 
are in the galaxy at the later epoch but are not in any re- 
solved galaxy at the earlier epoch as smoothly accreted gas. 
In Figure |3l galaxy mass is shown on the vertical axis, and 
time is shown on the horizontal axis. Each point represents 
an accretion event, showing the baryonic mass of the SKID 
identified galaxy when the particle was accreted onto it. Ev- 
ery particle that is in the z — group is shown. The locus 
traced by the topmost points shows the growth of the most 
massive progenitor of the z = galaxy, while lower ridges 
trace the growth of secondary progenitors which merge onto 
the main progenitor (or with each other) producing jumps 
in the mass track. There are occasional outlier points, cases 
in which an accreting particle is incorrectly assigned to a 
more massive neighboring system instead of the galaxy in 
which it is actually located. 

Following iKeres et al.) Hooi), we distinguish between 
cold and hot accretion. We trace the temperature his- 
tory of each particle and identify the maximum tempera- 
ture Tmax that the particle had prior to accreting onto a 
galaxy since the beginning of the simulation. Particles with 
Tmax < 2.5 X lO^K are classified as cold and shown in blue 
colour, and particles with Tmax > 2.5 x lO^K are classified as 
hot and shown in red. Since each point corresponds to accre- 
tion of one particle, intermediate mass galaxies sometimes 
show interleaved hot and cold accretion events. 

For satellite galaxies, we trace the history of its most 
massive progenitor galaxy and identify the epoch when it 
became a satellite, which we define as the epoch when the 
parent halo of the satellite merges with the parent halo of 
a more massive galaxy. In Figure |31 the magenta coloured 
vertical line marks the epoch when the galaxy became a 
satellite. 

Rows (a) and (b) of Figure [3] show high mass galax- 
ies with final baryonic masses (consisting of stars and cold, 
dense gas) of Mg ~ 2 — 3 x IO^'^Mq. Subsequent rows show 
galaxies with final baryonic masses Mg ~ 10^^ M© (c and d), 
Mg ~ 4 - 5 X 1O^°M0 (e and f) and Mg ~ 1 - 2 x 1O^°M0 (g 
and h). In each row, the first column shows a central galaxy 
in a halo of mass ~ 5 — 10 Mg. Subsequent columns show 
satellite galaxies in decreasing order of halo mass. The sec- 
ond column shows satellite galaxies selected from the most 
massive halo in the simulation {Mh ~ 3.4 x 10^** M©). The 
third and fourth columns show satellite galaxies from inter- 
mediate mass haloes and low mass haloes respectively. While 
these thirty-two examples cannot capture the full range of 
assembly histories that we see in the simulation, they suffice 
to illustrate the main qualitative trends. 

An examination of Figure[3]reveals a number of interest- 
ing trends. The assembly of high mass galaxies is complex, 
with active merger histories and continuing accretion at late 
times. Among high mass galaxies, there is not much differ- 
ence between the assembly histories of central galaxies and 
satellite galaxies, except for the satellites in the most mas- 
sive halo. Satellites in intermediate mass haloes, as shown in 



columns 3 and 4, continue to accrete gas and form stars after 
Zsat, the epoch when they become satellites. Some systems, 
such as the galaxy in the third column of row (b) and the 
fourth column of rows (a) and (c), "receive" mergers after 
becoming satellites. However, satellites in the highest mass 
halo are an exception. For these, assembly occurs early and 
there is little accretion after Zsat, and there is mass loss in 
some cases such as row (c) column 2. 

The assembly of lower mass galaxies appears to be 
smoother with fewer merger events. However, because our 
galaxy population is incomplete below our resolution thresh- 
old of 6.8 X 10^ M0, it is likely that we miss many lower mass 
merger events. Lower mass galaxies in low and intermediate 
mass haloes often have continuing accretion and star forma- 
tion after becoming a satellite, with row (g), column 3 being 
a striking example. Some lower mass galaxies in intermedi- 
ate mass haloes, such as column 3 of rows (e) and (f), have 
low but non-zero accretion rates at late times. The assem- 
bly histories of low mass galaxies in the most massive halo 
are systematically different from their central counterparts. 
They experience very little accretion after becoming satel- 
lites, and in many cases they experience sigiuficant mass loss 
(column 2 of rows f-h). 

Overall, Figure|3]shows surprisingly little distinction be- 
tween central and satellite galaxies. Galaxy mass is generally 
a stronger predictor of growth history than central/satellite 
status. We do see some trends of the expected sign, but Zsat, 
when the galaxy becomes a satellite, is not a sharp bound- 
ary in the accretion or merger history. Satellites in the most 
massive halo are an exception, in that accretion is slowed or 
shut-off at Zsa.t , and some satellites experience mass loss. We 
have examined assembly histories of many more systems in 
addition to the 32 shown in Figure |3l and they confirm the 
qualitative points made here. 

It is remarkable that satellites in low and intermediate 
mass haloes, which are moving at high velocities in the po- 
tential of the halo, continue to accrete gas and form stars. 
Our conjecture is that continuing accretion by satellites owes 
to substructure in the haloes. In effect, a satellite galaxy in 
a large halo may still be the central galaxy of its own co- 
herent sub-group. Figure |4] investigates one such system, a 
3 X 1O"M0 halo as it evolves from z = 1.07 to z = 0. 
We plot the dark matter particles in black and the baryonic 
particles in different colours, which depend on their galaxy 
membership at z = 0. Figure [S] shows the corresponding as- 
sembly histories of the eight most massive galaxies in this 
halo. Consider the most massive satellite galaxy, shown in 
red in Figure [3] (and in panel b of Figure [S| ■ It starts off 
as a central galaxy in its own halo. Following the merger of 
its parent halo with the parent halo of the central galaxy 
at z — 0.514, it continues to accrete gas from the region 
around it and has other galaxies gravitationally bound to 
it, in effect its own satellites. Even at 2; = 0, it is the cen- 
tral object of a distinctly defined dark matter substructure. 
Hence, accretion is not shut off at Zsa.t , and there is no sharp 
transition in its behaviour at this redshift. Another satellite 
galaxy, shown in green in Figure |3J becomes a satellite at 
z = 0.736, but it receives a merger after that time. 

We now turn to statistical measures to quantify the 
anecdotal results seen in Figures 3-5. Figure [S] compares the 
assembly histories of central and satellite galaxies as a func- 
tion of galaxy mass and parent halo mass. For four galaxy 
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Figure 5. Galaxy mass vs time for galaxies in the halo of Fig- 
ure |4] Panel (a) is the central galaxy, (fe) is the most massive 
satellite, shown in red in Figure l4l and (e) is the satellite shown 
in green in Figure |4] The format is the same as figure l3l 
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Figure 6. Fraction of final mass accreted by redshift z is plotted 
against redshift. We include mass in all progenitors of the 2: = 
system. In each panel, crosses show accretion histories of central 
galaxies, circles show all satellites, and lines show satellites in bins 
of halo mass as indicated in the legend of panel (a). 



mass bins, we plot the mass accreted by redshift z onto any 
of a galaxy's SKID-identified progenitors divided by the to- 
tal mass accreted up to 2 = 0. In each panel, the crosses 
show the results for all the central galaxies, the circles show 
all the satellite galaxies, and the lines show the satellites in 
different bins of halo mass. 

Looking first at central galaxies, there is a clear trend 



for higher mass galaxies to assemble at higher redshifts. For 

example, the redshift at which 50% of the final galaxy mass 

has been accreted by the central galaxy population (/acc=0.5 

in Figure [111 falls from z = 2.2 for = 10"-^Mo-10^^M, 
-,10.5 



2.2 for Mg = 1O"-^M0-1O 
to 2 = 1.2 for Ms = 1O^" ''Mq-1O^^M0, with a steady trend 







in between these masses. This downward shift of assem- 
bly redshifts is consistent with the well known "downsiz- 
ing" trend of observed galaxy evolution. It arises naturally: 
massive galaxies reside in dense environm ents with accel- 
erated early formation l|Neistein et al.|[2006l ). However, even 
the most massive galaxies have a significant fraction of mass 
accreted at low redshift, which is inconsistent with the ob- 
served red colours of the most massive galaxies. 

Low mass satellite galaxies accrete most of their mass 
early on and experience little accretion at late times, with 
/acc=0.5 at 2 = 2 compared to 2 = 1.3 for low mass centrals. 
There is a strong trend with halo mass: low mass satellites 
in the most massive halo assemble earlier and have accreted 
90% of their final mass by 2 = 1.5, while the satellites in 
Afh < 10^^ Mq haloes have slightly lower assembly redshifts 
than comparable central galaxies. The same trends are seen 
for higher galaxy masses, but the differences between central 
and satellite galaxies are smaller. For high galaxy masses 
(panels c and d), there are no satellites in low mass haloes. 

With our star formation prescription, which results in 
an increasing star formation rate with increasing gas density, 
star for mation closely trac ks the accretion of gas with a time 
lag (see iKeres et al.ll2008l ). We find similar trends for star 
formation as we did for gas accretion: satellite galaxies typ- 
ically form their stars earlier than central galaxies, but the 
difference is smaller at higher galaxy masses. Among satel- 
lites, those in high mass haloes typically form stars early on 
and have little star formation at late times, while satellites 
in low mass haloes continue to form stars at low redshifts. 

To further investigate the physical processes involved 
in gas accretion by satellites, we diffe rentiate between hot 
and cold accretion (|Keres et al.ll200^ ). There are two rea- 
sons that it is interesting to isolate the contr ibution of cold 
accr etion to galaxy gr owth. First, tests by iKeres I (|2007l ) 
and iKeres et al.l l|2008l ) show that the rate of hot gas accre- 
tion in hydrodynamic cosmological simulations is sensitive 
to differences among simulation codes but the cold accretion 
rates are robust. In particular, hot accretion rates a re lower 
in simulations with GADGET (jSpringel et aLlbOOll ) than in 
simulations with the PTreeSPH code used here, probably be- 
cause the former employs an entropy conserving formulation 
of SPH. Second, AGN feedback, which is not incorporated 
in our simulation, could be more effective at suppressing hot 
accretion flows, since they have low-densities and are quasi- 
spherical compared to col d accretion flows, which channel 
dense gas along filaments (|Keres et al.ll2005l ). Thus, for ei- 
ther numerical or physical reasons, the cold accretion rate in 
our simulation might be a better approximation to a galaxy's 
true total accretion rate. 

Figure [7] shows the fraction of 2 = galaxy mass that 
was originally acquired (by the galax y or its prog e nitors ) 
via cold accretion. In agreement with iKeres et al.l l|2005h . 
we find that low mass galaxies gain their mass almost en- 
tirely from cold accretion, while the highest mass galaxies 
gain up to 2/3 of their flnal mass by hot accretion. (This 
latter behaviour is different in the Gadget-2 simulations of 
iKeres et al.l[2008l ]. but the former remains true; see their Fig- 
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Figure 7. The fraction of 2 = galaxy mass that was gained 
through cold accretion, as a function of galaxy mass. The solid line 
shows central galaxies, the dotted line shows all satellite galaxies 
and other lines show satellite galaxies in bins of halo mass as 
indicated in the legend. 



Figure 8. Fraction of mass accreted in the cold mode vs redshift. 
This figure is analogous to Figure [6] but we eliminate all material 
accumulated through hot accretion. 



ure 10.) Figure [7] shows that the trend with galaxy mass is 
nearly identical for central and satellite galaxies, and that 
among satellite galaxies there is no dependence of the cold 
accretion fraction on halo mass. If hot accretion were com- 
pletely suppressed in the real universe (or in numerical simu- 
lations with a different treatment of gas physics or feedback) , 
then to a first approximation the z = galaxy masses would 
be reduced, on average, by the factor shown in Figure [T] 

Figure |8] repeats the assembly history analysis of Fig- 
ure |6] but we now eliminate all material gained originally 
through hot accretion (onto the galaxy or any of its progen- 
itors) before computing galaxy masses at any redshift. All of 
the trends seen in Figure|6]are still apparent: assembly shifts 
towards higher redshifts for higher galaxy masses, satellite 
galaxies assemble earlier than central galaxies of the same 
mass, and satellites in massive haloes assemble earlier than 
satellites in low mass haloes. However, the trend of assem- 
bly history with galaxy mass is much stronger (for centrals 
and satellites) because galaxies experience littl e cold accre- 
tion o nce their mass exceeds Mg ~ 3 x lO^^M© (|Keres et al.l 
I2OO5I ). Galaxies above Mg — 10^^ M© have typically accreted 
90% of their final mass by z = 2. Suppressing hot accretion 
is thus a way of turn i ng the most massiv e galaxies "red and 
dead" (lBin"nevll2004l: ICroton et al.ll200^: iKeres eralll2005l : 
iDekel fc Birnboimll2006l : ICattaneo et allbOOTf T 

So far, we have compared the properties of central 
galaxies to the properties of satellite galaxies over the course 
of their lifetime, from the time of their formation until z = 0. 
However, these satellites were once central galaxies, and they 
only become satellites at the redshift Zsat, when their parent 
halo merged with the parent halo of a more massive galaxy. 
To fairly compare the post-^aat growth of satellites to cen- 
tral galaxies, we create a "control" sample. We associate 
each satellite galaxy above the resolution threshold with a 
randomly chosen central galaxy of similar baryonic mass. 
We follow the evolution of this matched central galaxy from 
Zsat of the corresponding satellite galaxy to 2 = 0. Thus, 



the control sample has the same mass distribution as the 
satellites and is traced over the same range of redshifts. 

Figure |9] compares the mass accreted by satellite galax- 
ies after they become satellites to the mass accreted by 
the control sample of central galaxies over the same time 
period. Panel (a) shows the mass accreted by each satel- 
lite galaxy from Zaat to z = as a fraction of its mass 
at Zsat- Satellite galaxies in four halo mass bins are rep- 
resented by different colours and point shapes. For compar- 
ison, we show the mass accreted by our control sample of 
central galaxies over the same time interval. Central galax- 
ies tend to gain mass after Zsat of their matched satellites, 
with some even tripling their mass between Zaat and 2 = 0. 
For satellites, there is a trend with halo mass. Satellites in 
low mass haloes (M^ < 10^^ M©) gain a small fraction of 
their mass after Zsat- Satellites in intermediate mass haloes 
{Mh ~ 10^^ — tend neither to lose nor gain much 

mass; their mass at z = is similar to their mass at Zsat- 
Satellites in the high mass, Mh = 4 x 10^'' M© halo typically 
lose mass, in some cases up to half the mass they had at 
Zsat. The continuous nature of the trend with halo mass is 
more evident in panel (b) , which presents histograms of the 
accreted mass fractions summed over all galaxy mass bins. 

Panel (c) shows the mass accreted by galaxies in the 
cold mode from Zsat to 2 — 0, again scaled to the mass 
at Zsat (which still includes hot mode accretion up to that 
epoch). As in panel (a), we show satellite galaxies differen- 
tiated by halo mass, and we use the same control sample of 
central galaxies. For low mass galaxies, most of the accretion 
takes place in the cold mode. Hence, for low mass galaxies, 
panel (c) is identical to panel (a) . However, for higher mass 
galaxies most of the accretion takes place through the hot 
mode, so above Mg ^ IO^^'^Mq the accretion rates are lower 
for both central and satellite galaxies. As before, panel (d) 
presents histograms summed over bins of galaxy mass. It is 
again clear that central galaxies tend to gain mass after the 
matched Zsat and that satellites in the most massive halo 
tend to lose mass. However, with the lowered satellite accre- 
tion rates, the continuous trend in intermediate mass haloes 
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Figure 9. (Left) Mass change from accretion after becoming a 
satellite as a function of galaxy mass. Each point represents a 
galaxy. The red triangles, blue squares, green open circles and 
magenta filled circles show satellites in haloes of different masses. 
Each satellite galaxy is paired with a central galaxy of similar 
mass which is shown as a cross. Mg is the mass of the galaxy at 
2 = 0. Mzsat is the mass of all the progenitors of the galaxy at 
the epoch when the most massive progenitor became a satellite. 
Panel (b) shows a histogram of the fractional mass change from 
accretion for satellites in different halo mass bins as well as central 
galaxies. (To preserve visual clarity, we plot this histogram as 
connected points.) (Right) Mass change from cold accretion after 
becoming a satellite. Panel (d) shows a histogram of the fractional 
mass change from cold accretion for satellites in different halo 
mass bins as well as central galaxies. 



is no longer so clear - there is only a small amount of satel- 
lite accretion in the cold mode, even in relatively low mass 
haloes. 

In our preceding analysis, we examined the mass ac- 
creted from Zsat to the present time and found differences in 
the total amount of accretion between central and satellite 
galaxies. To examine the timescale over which these differ- 
ences develop, we define a dimensionless accretion rate as 
follows: 



n = 



AMacc H' 



At 



(1) 



where AMacc is the mass accreted in time interval At, H is 
the Hubble parameter evaluated at the centre of the time 
interval, and Mg is the mass of the galaxy at z = 0. For 
each epoch, we take a small time interval around the cen- 
tral epoch and calculate the accretion rate using the mass 
accreted, AMacc and the time interval At between those 
epochs. The precise interval At depends on the simulation 
outputs available near Zaat, but typically At ~ 300 Myr. 

Figure [To] compares the accretion rate of satellite galax- 
ies to the accretion rate of central galaxies in the control 
sample at four different epochs. Panel (a) shows the accre- 
tion rate at Zsat when the galaxy becomes a satellite, while 
panels (b), (c) and (d) show the accretion rates 200 Myr 
after z^at, 500 Myr after Zsat and f Gyr after Zsat, respec- 
tively. The median Zsat values in the four halo mass bins are 
0.19, 0.55, 0.59 and 0.99 in increasing order of halo mass. 
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Figure 10. Accretion rate against galaxy mass in bins of 0.25 
dex. The top left panel shows the accretion rate at Zsat , the epoch 
at which the galaxy becomes a satellite. The other panels show the 
accretion rate 200 Myr, 500 Myr and 1 Gyr after Zsat. Accretion 
rate is defined as, rj = (AMacc/Mg) (H~^ /t) where Mace is the 
mass accreted in time t, is evaluated at the centre of the time 
interval, and Mg is the mass of the galaxy at z = 0. Lines show 
the average value for all satellites (dotted) and matched central 
galaxies (solid) in each galaxy mass bin. Points with error bars 
show the average value and uncertainty in the mean in four bins 
of halo mass, as indicated in the legend. 
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Figure 11. Cold accretion rate vs galaxy mass. This figure is 
analogous to Figure [TOl but showing only cold accretion. 
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The dotted and solid curves, averaged over all satellites and 
matched centrals in each galaxy mass bin, afford the best 
statistics for comparison. At Zsat and Zsat+200 Myr, the di- 
mensionless accretion rates of satellite and central galaxies 
are similar. After 500 Myr, the accretion rates of satellites 
are systematically lower, and the difference increases at 1 
Gyr. Only the lowest mass satellites have significantly non- 
zero accretion rates at z^ax + 1 Gyr. Points with statistical 
error bars show the accretion rates in bins of halo mass. The 
drop in satellite accretion rates develops most clearly in the 
highest mass halo, followed by the next highest mass bin. 
By Zsat + 1 Gyr, some galaxies in the highest mass halo are 
losing mass. While this analysis would certainly benefit from 
better statistics (i.e. a larger simulation volume), the trends 
with time and halo mass are continuous. They indicate that 
satellite accretion shuts off gradually, over a timescale of 0.5- 

1 Gyr, with low mass galaxies in high mass haloes showing 
the strongest effect. 

Figure [TT] repeats this analysis for cold accretion only. 
Here the accretion rates are close to zero for all galaxies with 
Mg above ~ lO^^'^M©, both central and satellite, at all the 
epochs shown. At lower Mg, the accretion rates are still sig- 
nificantly reduced relative to Figure [TUl but the differences 
between central and satellite galaxies and the dependence 
on halo mass are similar. It is particularly notable that low 
mass satellites in the IO'^^-IO'^^Mq halo mass bin have pos- 
itive cold accretion rates at Zsat + 200 Myr and Zsat + 500 
Myr. Note, however, that the most massive haloes in the bin 
are ~ 3 X IO^^Mq. Even satellites in the 3 x lO^^M© halo 
show cold accretion at Zsat and, to a lesser extent, at Zsat + 
200 Myr. 

Besides smooth accretion of gas, galaxies also grow by 
receiving mergers from lower mass galaxies. Since we do not 
resolve galaxies below 6 x lO^M©, we restrict our merger 
analysis to systems in which the primary galaxy is more 
massive than 3 x 10^° and the merger mass ratio is 
i?mcrg > 0.25. As we go to higher redshifts, the galaxy 
masses decrease, causing the number of resolved mergers to 
decrease. Hence, we restrict our merger analysis to mergers 
between z = 1 and z = 0. The top panel of Figure [T^ shows 
JJmerg for all merger events between z = 1 and 2: = as a 
function of the z = G baryonic mass of the galaxy, differen- 
tiating between the central/satellite status of the primary 
galaxy at Zmorg (the epoch of the merger) and at 2 = 0. We 
see mergers with a wide range of mass ratios for satellites 
as well as centrals. A small fraction of mergers (~10%) are 
between objects of similar mass (i?morg > 0.9). More than 
four-fifths of merger events involve a satellite galaxy merging 
with the central galaxy of its parent halo. About one-third 
of these galaxies go on to become satellites by z = 0, follow- 
ing a merger of their parent halo with a halo hosting a more 
massive galaxy. However, it is notable that close to one-fifth 
of the merger events in our simulation are mergers between 
two satellite galaxies. 

The bottom panel of Figure [12] compares the fraction of 
central and satellite galaxies that receive mergers between 

2 = 1 and 2 = above a certain merger ratio, Rmag. One- 
fourth of the galaxies in the simulation (above 3 x 10^" M0) 
experience at least one merger with Rmerg > 0.25 between 
2 = 1 and 2 = 0. This fraction is similar for central galaxies 
and for 2 = satellites and the Rmag distribution is also 
similar, as one can see by comparing the solid and dashed 
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Figure 12. (Top) Merger events after 2 = 1 with mass ratio 
RmcTg > 0.25, versus galaxy baryonic mass at 2 = 0. Central 
galaxies, 2 = satellites that are centrals at the merger epoch 
Zmcrg and galaxies that are satellites at Zmerg are shown sepa- 
rately, as indicated in the legend. (Bottom) Fraction of central 
galaxies (solid line), 2 = satellites (dashed line) and galaxies 
that are satellites at 2incrg (dotted line) that receive a merger 
above a certain -Rmorg- 



lines. However, the fraction of satellites that receive merg- 
ers after 2sat is a factor of 2-3 lower (dotted line). Galaxies 
that receive mergers often receive more than one merger, so 
the accounting in the top and bottom panels is somewhat 
different. 

Figure [12] adopts the somewhat arbitrary redshift inter- 
val 2 = 1 to 2 = 0, and the mass distributions of the central 
and satellite galaxies are systematically different. To fairly 
compare the merger activity of central galaxies and post-2sat 
satellites, we return to the control sample used for our ac- 
cretion comparisons. The top panel of Figure [13] shows the 
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Figure 13. (Top) Merger mass ratio -Rmcrg vs 2 = galaxy 
baryonic mass for mergers received by satellite galaxies after Zsat 
compared to the control sample of central galaxies over the same 
time period. (Bottom) Fraction of satellite galaxies (dashed line) 
and central galaxies in the control sample (solid line) that receive 
a merger above a certain -Rmcrg- 



merger ratio for mergers (no longer restricted to z < 1) with 
J?morg > 0.25 received by satellite galaxies and by the control 
sample of central galaxies after Zsat of the matched satellite. 
Only 14% of satellite galaxies and 16% of central galaxies 
in the control sample receive mergers after Zsat- The dis- 
tribution of mass ratios is broad for both populations. Our 
sample is not large enough to study the merger properties 
of satellite galaxies as a function of halo mass. The bottom 
panel of Figure [13] shows the fraction of central and satel- 
lite galaxies that receive a merger above a given 
Zsat (in contrast to Figure W2l which shows all mergers since 
2 = 1 and includes all central galaxies) . The probability of a 
satellite galaxy receiving a merger above a certain mass ra- 



tio is only ~2% lower than that of a central galaxy of similar 
mass, across a wide range of merger mass ratios. 



4 DISCUSSION 

The central galaxies of haloes form a distinct population 
in t his simulation, as an t icipated by semi-analytic mod- 
els llWhite fc FrenkI Il99ll : iKauflFmann. White fc Guideronil 
ll99a : ICole et al.lll9 94'~ and as found in previous numerical 
studies (jBerlind et al . 2003: Zheng et al. 200_5). In the great 
majority of FOF haloes, the most massive galaxy lies at the 
bottom of the dark matter potential well and close to the 
centre of mass. In more than two-thirds of the FOF haloes, 
the baryonic mass (stars -I- cold gas) of the central galaxy 
is 2-5 times higher than that of the most massive satellite. 
Only a small fraction of FOF haloes host two galaxies of 
similar mass. Because of its limited volume, our simulation 
has only one cluster mass halo {Al ~ 3 x 10^*Mq). The 
other haloes massive enough to host resolved satellites range 
"M© to ~ 3 X 10^^ M0, while haloes down to 



from ~ 3 X 10 

~ IO^^A/q host resolved central galaxies. 

The distinction between central and satellite galaxies in 
our simulation is weaker than expected in a simple picture 
where only central galaxies accrete mass and "receive" merg- 
ers of less massive systems. Galaxies that are satellites at 
z — assembled their mass earlier, on average, than present 
day central galaxies of the same baryonic mass, and satellites 
of more massive haloes have systematically higher assembly 
redshifts. However, to a large extent the baryonic mass of a 
galaxy is a better predictor of its accretion and merger his- 
tory than its central or satellite status, with the exception 
of galaxies in the cluster mass halo. In part, the dominance 
of baryonic mass over satellite status just reflects the fact 
that many satellites were central galaxies for most of their 
history, merging into larger haloes only at low redshift. How- 
ever, we find that simulated galaxies continue to accrete gas 
and, in some cases, merge with lower mass objects even af- 
ter they become satellites. In essence, galaxies that enter 
the virial radius of a larger halo may remain the "central" 
galaxies of their own dark matter substructure, and they lose 
their central status only gradually. Alternative halo defini- 
tions (e.g., higher threshold density or spherical overdensity 
identification) could yield somewhat sharper distinctions be- 
tween central and satellite galaxies, but we do not expect 
them to change this basic result. 

While galaxies continue to accrete gas after becoming 
satellites, they accrete less than a "control" sample of cen- 
tral galaxies with similar baryonic masses measured over 
the same time period. At the time they enter a larger halo, 
satellites have similar accretion rates to central galaxies of 
the same mass, and the accretion rate declines towards zero 
over the next 0.5-1 Gyr. Lower mass satellites in the cluster 
mass halo begin to lose mass, and by z = they are typically 
10-20% less massive than they were when they became satel- 
lites. When we compare satellite galaxies to central galaxies 
of the same baryonic mass, we flnd large systematic differ- 
ences if, and only if, the mass of the satellite's host halo is 
much larger than the typical halo mass for the corresponding 
central gal axies. 

As in iKeres et al. 

I (12005'), we flnd that galaxies below 
""Mq accrete predominantly through "cold 
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mode" , while higher mass galaxies grow mainly through "hot 
mode" accretion of virial temperature gas. This transition 
is similar for central and satellite galaxies. Consequently, 
our general conclusions would still hold if we assumed that 
a galaxy's true accretion rate was close to its cold (rather 
than total) accretion rate, for either the physical or numer- 
ical reasons discussed in §3, except that the accretion rates 
of massive galaxies would become very low. Regardless of 
whether we consider total accretion or cold accretion only, 
these simulations exhibit natural "downsizing" in which ac- 
cretion and sta r formation shift toward s lower mass galaxies 
at later times (|Keres et al.ll2005l . bOOSj ). 

Comparing present day central galaxies and present day 
satellites, above a mass threshold of Mgai = 3 x 10^" Mq, we 
find that only a slightly larger fraction of centrals (27% vs. 
22%) have received mergers with > 0.25 

since z = 1. Roughly half of the satellite mergers occurred 
after the larger parent had become a satellite, and half oc- 
curred while it was still a central object. The mass ratio 
distributions are similar in all cases. There are more central 
galaxies than satellites and only about one in five resolved 
merger events involves two satellites. However, if we com- 
pare satellites to a control sample of centrals as we did for 
accretion, we find similar fractions experiencing resolved. 



> 0.25 mergers during the course of the simulation. 



with similar Rmerg distributions. In summary, most galaxies 
have quiescent assembly histories at z < 1, and the merger 
rates of central and satellite galaxies are not radically dif- 
ferent. The most massive galaxies do have more complex 
merger histories (|Maller et al.ll2006l ). and these are usually 
central galaxies of high mass haloes. 

While we have focused on mass growth rather than 
star formation, the two are closely linked in our simula- 
tion, with star formation generally following gas accretion 
after a short delay. The fact that galaxies experience con- 
tinuing gas accretion for 0.5-1 Gyr after becoming satellites 
means that satellites will be systematically bluer than pre- 
dicted in semi-analytic models that assume no accretion onto 
satellites ( "strangulation" ) . There are numerous indications 
that including the continuing satellite accretion predicted 
here would improve t he match between semi- analytic mod- 
els and observation s. IWeinmann et all (j2006h find that the 
ICroton et aP (2006) semi-analytic model predicts an exces- 
sive fraction of red satellites relative to the satellite popula- 
tion in groups identified from the Sloan Digital Sky Survey 
(SDSS). The discrepancy is largest for lower mass groups, 
which have much higher blue fractions than predicted, con- 
sistent with our finding that the central-satellite di stinction 
becomes sharper for satellites in high mass haloes. iLi et al.l 
l|2007l ) show t hat the pai r wise veloci ty distributio n s pre - 
dicted by the iKang et al.l (|2005| ) and ICroton et all (120061 ') 
models can be brought into agreement with SDSS mea- 
surements by re ducing the fraction of faint red galaxies. 
ICoil et al.l ll2008t) find s imilar trends at 2; 1, showing that 
the lCroton et al.l l|2006l ) model predicts excessive clustering 
of red galaxies and insufficient clustering of blue galaxies 
relative to the DEEP2 Galax y Redshift Survey, esp ecially 
on small scales. Most recently, Weinmann et al.l l|2008l ) show 
that the colours of SDSS satellites can be reproduced if their 
star formation slows over a timescale of 2-3 Gyr after Zsat, 
somewhat larger than the timescale for accretion shut ofi' 
found here. 



Other observations indicate that the colours of satel- 
lite galaxies depend systematically on halo mass, as we 
would expect based on the assembly redshift trends in Fig- 
ures 6 and 8. In a sample of clusters tak e n from the Red- 
Sequence Cluster Survey, ICilbank et al.l (|2008l ) find that 
the red fraction in clusters increases with richness, indicat- 
ing an earlier trunca t ion o f star formation in higher mass 
haloes. [Hansen et all (|2007l ) find a similar result using data 
from the SDSS. Finally, X-ray studies of clusters indicate 
that massive satellites retain local concentrations of h ot gas 
(jSun et all I2OO7I : [ Jeltema. Binder fc Mulchaevl[2008f ). pro- 
viding further evidence that the physical transition from 
"central galaxy" to "satellite" is a gradual one. Further dis- 
cussions of the impact of the halo environment on satel- 
lite galajcy properties, focuse d on the effects o f ram pres- 

Dpin 

iFont et al.l l|2008l ). 

Our small simulation volume limits our statistical 
power, especially for massive haloes. In addition, there are 
indications that the TreeSPH code used here predicts too 
much hot gas accretion, and that the entropy conserving 



sure stripping, a re given by iMcCarthv et al.l (2008) and 



SPH formulation o flSpringel fc Hernguistl (12002 ) yields more 
accurate results (ISpringel fc Hernguistl I2OO2I : iKereF et al.l 
l2008l : lKeres1l2007l ). Keres et all (|2008D have recently anal- 
ysed a simulation of a 50/i~^Mpc cube ev olved with the 
entro py conserving SPH code GADGET-2 (|Springel et all 
l200lD . with a similar mass resolution to the simulation anal- 
ysed here. While they focus mainly on the relative impor- 
tance of cold and hot accretion in galaxy growth, they also 
investigate central and satellite accretion, and confirm our 
main result: satellites experience significant amounts of con- 
tinuing accretion, and the distinction between central and 
satellite growth rates is smaller at high galaxy masses and 
(a point not invest igated here) at hig h redshift. Both our 
simulation and the iKeres et al. I (|2008') simulation overpre- 
dict the observed galaxy baryonic mass function, probably 
because they do not include galactic winds that eject ac- 
creted material from low and intermediate mass galaxies. 
Winds will change the relation between accretion rates and 
star formation rates, and they will affect the accretion rates 
themselves by altering the local gas supply. However, we 
do not expect them to alter our qualitative conclusions: the 
distinction between central and satellite galaxies is funda- 
mental to understanding the galaxy population and its evo- 
lution, but differences in their growth rates and hence their 
physical properties are smooth, not sharp. 
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